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A family of multiple integrals over four variables is rewritten in terms of a family 
of simple integrals involving the product of four modified Bessel (Macdonald 
functions). The latter are shown to be related to 7£(3)/2. A generalization to 2n 
, integration variables is given which yields only £ at odd arguments. 
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abc + bed + cda + dab 
where a = fa — fa, b = fa — fa, c = fa — fa, d = /3 + /3 4 — fa. One has 

r /9-a ^-«-6 rf (a + c)(& + d) 
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hMP)= da db dc- 
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^ . Jo Jo Jo a bc + &cd + cda + da& 
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where d = (3 — (a + b + c) is understood 1 . By successive integration on c, b, a 
one has obtained [I] 

Id 

where C(3) = 7C(3)/2. Rewrite Jnoo(/9) differently 

W0) = I" dadfrdedd - ^t^ 6 /?, , ^-(a + 6 + c + d)) 
Ja,b,c,d=o a bc + bed + cda + dab 

«/?)=/ dadfcdc + ^ + 

JO<a+f)+c<0 



abc + bed + cda + dab 



and also 

hwo(P) =(3 dadbdc dd 5({3 -(a + b + c + d)) 

Laplace transform 

1100 ^ ^ Jo ^ P s 3 J a ,b,c,d=o a ^ c + bed + cda + dab 

such that /noo( s ) = ^hioo(P = l)/ s3 - Fix s — 1 and denote /noo( s — 1) — 
Jiioo- One has 

/*oo POO 1 

7noo= / dadfedcdd / dt^e^ a+b+c+d ^ + ^+^ 

Ja,b,c,d=0 JO C d 

Change variable u = 

where the K^uYs are modified Bessel functions also called Macdonald func- 
tions [2 . Therefore 

/no. = 2< [ ^(IWW = 2 i±M (1) 
2 Generalization 

To generalize (JTJ), introduce for n a , rib, n c , rid E N 

r°° n n a hn b n c Jn d 

/nn„nn„= / da db dc dd - , e -(a+&+c+<*) (2) 

Ja,b,c,d=0 a " C + " C " + C " a + « a " 

and rewrite it as 2 
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2 This is also formally true for n Q ,n(,,n c , £ Z but diverges as soon as n a < or 
Ub < or n c < or < 0. For n a , nb,n c , rid. > m — 1, to > 1 



(a6c + 6cd + cda + da&) r 



IZn b n c n d = / da d6 dc dd , , " ^ " I " - , , ; ,„ e ~^^ 



One has 
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n a n b n c n d + Vn a n b n c nM (3) 



- L n a n b n c n d 

where u na n bncnd , v nanbUcnd are positive or negative rational numbers. 
-Proof: Since 

K na (u) = 2(n a - l) Kn ^ u) + K na _ 2 (u) 

u 

then if n a > even/odd 



u 

+ CiuKiiu) 
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^K na (u) = Z-K Q/1 (u) + ^—- K 1/0 (u) + c na _ 2 u na 2 K 0/1 (u) + ... 



where c na _ 2 , . . . , c\ are rational numbers. It follows that © is true if 
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p n (0000)= / duu n+1 K {uf n>0 (4) 
Jo 

POO 

p n (0011)=/ duu^K^ufK^uf n>2 
Jo 

POO 

p n (llll)= / duu^K^uf n>A (5) 



in even 



rewrites as 



L ^ ( " r „+„ 6+ „ c+ „ d -4( m -l)^^^ 



r(m) J 2 v 2' 

The considerations that follow for I na n b n c n d (i.e. the case m = 1) generalizes to I™ a n b n e n d 
for m > 1 implying that the latter can also be built from p„(0000) and p n (llll) defined in 
iJU and (JSJ. Therefore I™ nbncni is of the form (J3J). Note that if one introduces back s in 

jm 

n a n b n c n d 

T m ( s )= / dad&dcdd - s (a+f,+c+d) 

l*J y O)6fCj|J=0 d ° M aC ad (abc + bed + cda + dab)™ 6 

and differentiates n times with respect to s, one obtains with s = 1 

r°° n n a f,n br n cr in d 

C»>n»* = / dadbdedrifa + fr + c + o')" ° C ° , e -< Q+b+c+d ) 

™ a " i> " <= " d J a ,b,c,d=o {abc + bed + cda + dab) m 

again of the form |J3}. 



and 

i„(0001) = / duu^K^ufK^u) n>\ 
Jo 

POO 

z n (0111)=/ duu n+1 K Q {u)K x {uf n>3 
Jo 

(n odd) 

are again 3 of the form (J3J). To show this, integrate by parts (use dK (u)/du = 
—Kxiu) and d{uKi{u)) / du = —uK (u)) to obtain 

p n (0000) = — ^-z n+1 (0001) n > 
n + 2 

p n (0011) = -(i n+ i(0001) + i n+ i(0111)) n > 2 
p n (llll) = — ^-z n+ i(0111) n > 4 

which implies 

2np n (0011) = (n + 2)p„(0000) + (n - 2)p„(llll) n > 4 
It remains to be found p n (0000) and p n (llll): integration by parts gives 

z n (0111) = — ^-((p„ +1 (llll) + 3p n+1 (0011)) n > 3 
n — 1 

, n (0001) = ^-(( Pn+1 (0000) + 3p n+1 (0011)) 7i > 1 

71 + 1 

thus the recurrence relation acting in a 2 dimensionnal vector space 

p n+2 (0000) \ _ 1 / (n + 2) 2 (5n + 4) -3n 2 (n - 2) \ / p n (0000) 



p„ +2 (llll) J 2 5 (n + 2) V -3(n + 2) 2 (n + 4) n(n - 2)(5n + 16) ^ V Pn( im ) 

(6) 

n > 4. 

-Examples: 



3 Note that a product of four Macdonald functions is needed here: for example 

poo Q 00 1 o 



p=0 

is not simply given in terms of £'s. 



by direct computation < n < 4 

C(3) 



Po(0000) 
i x (0001) 
p 2 (0000) 



2 2 

cm 

2 3 

-3 + C(3) 
2 4 



p 2 (0011) = 1 + ^ known[U 

i3( oooi) = dLtM 

^(Olll) = — obvious 
and n = 4 (initial conditions) 

-3 3 + 7C(3) 



p 4 (0000) 

P4(1H1) 



2 6 

53 - 3 2 C(3) 
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then the recurrence (JBJ) gives for n > 4 

-37 + 3 2 C(3) 
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2 4 

, 67-(3)(5)C(3) 

2^ 

-(5)(19)(269) + (3 2 )(7)(97)C(3) 



2 io 

nn) = 3 (13)(811) -(3 2 )(5 2 )(11)C(3) 



Pio(OOOO) 
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9304913 + (3 2 )(5 3 )(11)(179)C(3) 

WW) 



„ ri111 x o4 (3)(ll)(4139)-(5 3 )(7)(37)C(3) 
p 12 (0000) = 3 



, (7)(19)(23909) + (5 3 )(23)(263)C(3) 



(2 12 )(5) 

„r43Vfi7Vl15l! 

p 12 (llll) = 3 



2 (43)(67)(115 19) - (3 2 )(5 3 )(7 2 )(11)(13)C(3) 



(2 12 )(5) 



If one sets n = 2k and defines q k (0) = p 2k (0000)/(2k)\, q k (l) = p 2fc (llll)/(2fc)! 
the recurrence (jUJ) becomes 

g fc+ i(0) \ _ 1 ( (5A: + 2) -3A; 2 (A;-1)/(A; + 1) 2 \ f q k (0) 
g fc +i(l) J 2 4 (2A; + 1) V ~3(A; + 2) fc(fc - l)(5fc + 8)/(fc + l) 2 J \ q k (l] 



k > 2. In the asymptotic regime k — > oo the recurrence matrix 



(7) 



1 ( 5 " 3 
2 5 V -3 5 



diagonalizes as 



with eigenvalues 1/4 and 1/16. In this convenient diagonal basis (J7|) becomes 

9fc+i(l) +9fc+i(0) 
g fe +i(l) - 9fc+i(0) 

1 / A; — 2 + yfc(fc — l)(fc + 4)/(fc + l) 2 -(Jfe - 2) + k{k - l)(k + A)/(k + if 



2\2k + 1) V 4 (-( fc + 1) + k(k - l)/(k + 1)) 4(fc + 1 + fc(fc - + 1)) 

gjb(l) + gjb(O) 
<7*(i)-g*(o) 

fc > 2 with by definition q k {\) + g fe (0) > g fc (l) - q k (0) > 0. 
-Examples: 
k = 2 (initial conditions) 

?2(l) + ?2(0) = ^y|y^ (8) 

*(U-*«» = ^f o» 



fc > 2 

g 3 (i) + ?3(o) 
«3(i)-? 3 (o) 



53 - (3) 2 C(3) 
(2 10 )(3 2 )(5) 

349 - (3) 4 C(3) 
(2 10 )(3 2 )(5) 





+ 94(0) = 


3037- (3) 2 (73)C(3) 
(2")(3*)(5)(7) 


94(1) 


-94(0) = 


1787- (3) 2 (47)C(3) 
(2 12 )(3 2 )(5)(7) 
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+ 9s(0) = 


(439)(2003) - (3 2 )(5 3 )(181)C(3) 
(2")(3*)(5»)(7) 


95(1) 


-9s(0) = 


(7)(73)(1993) - (3 2 )(5 4 )(43)C(3) 
(2^)(3 4 )(5 2 )(7) 
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+ 9e(0) = 


(2283583) - (3 2 )(5 3 )(479)C(3) 
(2 21 )(3 3 )(5 3 )(7)(11) 


9e(l) 


-9e(0) = 


(127)(1901)-(3 3 )(5 3 )(17)C(3) 
(2")(3»)(53)(7)(11) 


9r(l) 


+ 97(0) = 


(53) (1708543) - (3)(5 3 )(7 3 )(167)C(3) 
(2 2 i)(3 2 )(5 3 )(7 3 )(ll)(13) 


97(1) 


-97(0) = 


(13) (61485173) - (3)(5 3 )(7 4 )(211)C(3) 
(2 23 )(3 2 )(5 3 )(7 2 )(11)(13) 



Both g fc (l)+g fc (0) and g fc (l)-g fe (0) decrease (faster than (l/16) fe ~ 2 ) as (l/16) fc_2 
(when k finite) when k — > oo, i.e. their asymptotic behavior is governed by 
the smallest of the eigenvalues. This is due to the initial conditions fine tun- 
ing (JHJ) and 0, which defines a £(3) dependant initial orientation in the 2 
dimensional vector space leading to this particular asymptotic behavior. This 
situation is analogous to the one encountered for establishing the irrationality 
of £(3) (see [3] for a review of various demonstrations of the irrationality of 
C(3)): in that case a second order recurrence relation (or equivalently a family 
of integrals) with £(3) dependant initial conditions such that an asymptotic 
convergence governed by the smallest of eigenvalues is achieved. It is how- 
ever manifest in the present case that both q k (l) + q k (0) and q k (l) — 9fc(0) 
have denominators which increase too fast with k to obtain a convergence al- 
lowing for an alternative proof of the irrationality of C(3). This is still the 
case for any fc-independant initial conditions rotation in the plane. If one 
allows for more general fc-dependant combinations, note that if one redefines 
9fe(01) =p2fc(0011)/(2fc)! then 



q k (01) = q k (l) + q k (0) - (q k (l) - q k (0))/k 



if one redefines j' fc (0111) = z 2fc+ i(0111)/(2A; + 1)! then 

and finally if one redefines j fc (0001) = i 2fc+ i(0001)/(2fc + 1)! then 

^< 0001 > = 2^FTi)*( 0) 

Particular &;-dependant linear combinations of ^(1) + ^fc(O) and gfc(l) — gfc(O) 
follow as 

A(oui) + i 4 (oooi) = < *(D + *«>) 

i»(oni) - A (oooi) = tWD-^W-y 1 ) + «*(")) < a( i) _ a(0) 

with by definition j fc (0111) - j fc (0001) > 0. However both j fc (0111) -^(0001) 
and jfc(Olll) + jfc(OOOl) have again denominators which increase too fast with 
k. 

3 Further Generalization: C(5), ((7), ■ • • 

In view of obtaining ((5), a natural generalization of (J2J) which would still 
involve Macdonald functions would be, in the simplest case 

/ da db dc dd de df — — — — — — — ^_^ e -(a + b + c + d + e+f ) 

J a ,b,c,d,ej=o abode + bcdej + cdeja + dejab + ejabc + jabcd 



nothing but 



2 6 / ^K (uf 



o 



Numerics indicates however that this integral is not a linear combination with 
rational coefficients of 1, C(5), C(3)> C(2)C(3)> • • • • Introduce rather for n > 2 
the 2n integration variables a, and 6j, with < % < n, and consider 



poo n n 

/ |T ddi dbi 5(1 - YVoj + bi)) — - 
Jo £5 £i ai6ii 



P(ai,bi, . . . ,a n ,b nJ 



t | , _ x (a 2 + b 2 ) ... (a n + b n ) + cp 



where cp means cyclic permutation over 1, 2, n and P(ai, 61, ... , a n , b n ) is 
any given monomial of the variables a i; 6j. The case n = 2 corresponds to the 
variables a, 6, c, d used up to now. This is a family of integrals with root 



poo "> ,b 
In= IT d °* db * <5 ( 1 ~ + T 



Change variables 
to get 



Ui = a* + bi Vi 



(a 2 + b 2 ) ... (a n + b n ) + cp 
as - bi 



a-i + k 



/»oo n n „i n ... 

7 " =4 y. nd».*(i-E».)y o rK (1 _ 1 ,, )tI1+ ... +(1 _ t , ; 



)Ur. 



In order to integrate over u\. . . u n , follow the procedure used above: Laplace 
transform the S measure and exponentiate the denominator 



In = 



r(n- 1) 

Integration over Ui and Vi leads to 



pi n poo n / n \ 

y JJd^ y dt JJd^ exp f -^^(l + t(l-^ 2 ))J 



1)/ 



dt 



which can be rewritten 4 as 

2™ 



(v / 0+^) n 



1 dy / j/ 



log 



n— 1 



+ 1 



1+t _ 1 
t ' 



Using 



leads to 



*or as 



r(n-i)7 y V 1 -!/ 2 

rfn + i 



(log-)" 

y 



fdyy^log-r 

Jo y 



(p + i) n+1 



In = 2 



n r(n+l) ~ (p+l)(p + 2)...(p + n-2) 
r (™ - !) 2 ^ (n + 2p - 



(10) 



In = 



n 



V sn-l 



-Examples: 

h = 7C(3) 
h = 3C(3) 



31 

A = 3(7C(3) - -C(5)) 



h = |(C(3) - C(5)) 

h = ^(TC(3)-^C(5) + ^C(T)) 

Only C(°dd) does appear in ffTTTft as can be seen by distinguishing between n 
odd and n even 

r _ 4r(n + 1) ^ rfc 2)/2 ((^ + 2p - l) 2 - (2g - l) 2 ) 
n ~ I>-1) 2 ^ (n + 2p-l)»+i ^ 6Ven 

4r(n + 1) ~ nfc 3)/2 ((n + 2 P - l) 2 - (2g) 2 ) 

i„ = — r— > ; : TL Odd 

T(n-l) 2 ^ (n + 2p-l) n 
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